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ABSTRACT 

We present the simultaneous Swift and Fermi observations of the bright GRB 100728 A and its 
afterglow. The early X-ray emission is dominated by a vigorous flaring activity continuing until 
Iks after the burst. In the same time interval high energy emission is significantly detected by 
the Fermi/LAT. Marginal evidence of GeV emission is observed up to later times. We discuss the 
broadband properties of this burst within both the internal and external shock scenarios, with a 
particular emphasis on the relation between X-ray flares, the GeV emission and a continued long- 
duration central engine activity as their power source. 
Subject headings: gamma-ray burst: individual (GRB 100728A) 
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1. INTRODUCTION 

The Fermi Gamma-Ray Space Telescope, launched in 
June 2008, has taken the study of GRBs into an energy 
realm that so far has b een poorly explored . Fermi /h AT 
(Large Area Telescope: I Atwood et al.|[2009D observations 
of GRBs allow for the first time a detailed study of 
the temporal and spectral behavior at high energies 
(>100 MeV). One of the most interesting features is 
the detection of a delayed and rapidly decaying high- 
energy emission, lasting hundreds to thousands of sec- 
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onds longer than the observed sub-MeV 7-ray emis- 
sion (|Abdo et al.l [2009a. b.c). Extended GeV emission, 
first hinted at in EGRET observations ()Hurlev et al.l 
[T99l . appears now as a common feature of Fermi /L AT 
bursts. The nature of such long-lived high-energy emis- 
sion is far from being established. One possibility is that 
it is generated via synchrotron radiation of the exter- 
nal forward shock (iKiimar fc Barniol Duranll2009l 120101: 
iGhisellini et al.l I2O1O0 . An alternative scenario is that 
it refiects the gradual turn-off of the central engine ac- 
tivity (Zhang et al. 2010). Such interpretations predict 
very different afterglow behaviors (jPiran fc Nakadl26l0l : 
iMimica et al.l [201011 and therefore can be directly veri- 
fied through broadi)and (from optical/X-ray to GeV en- 
ergies) early-ti me observations. To date , only one burst 
(GRB 090510: [PePasquale et ;iI1[2010l ) of the 20 LAT 
detected GRBs has been simultane ously detected by th e 
Swift multi- wavelength observatory (|Gehrels et al.l200^ . 
In this case, an afterglow emission pro vides a likely expla- 
nation of the broadban d dataset (e.g. iDe Pasauale et all 
[MotlCOTsTeFalllMol) . 

In this Letter, we report on the Fermi /h AT detection 
of a temporally extended emission from GRB 100728A 
and the simultaneous Swift observations of an intense 
X-ray flaring activity. We further discuss the possibil- 
ity that in the case of GRB 100728 A the observed high- 
energy emission is related to X-ray flares and ultimately 
to the long-lasting activity of the inner engine. Obser- 
vations and analysis are reported in § [21 our results are 
discussed in § [31 we draw our conclusions in § [H Unless 
otherwise stated, the quoted errors are at the 90% con- 
fidence level and times refer to the Fermi /GBM trigger 
To. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Swift data 

The bright GRB 100728A came into the Swift field 
of view during a slew to a pre-planned target, when 
the trigger system is disabled. After the spacecraft set- 
tled, the burst triggere d the Burst Alert Telescope (BAT; 
iBarthelmv et al.ll2005[ ) on-board Swift at 02:18:24 UT on 
28th July 2010. Swift slewed immediately to the burst. 
The t wo narrow field instr uments, the X-ray Telescope 
(XRT; IBurrows et aI[|2005D and the Ultraviolet Optical 
Telescope fUVOT: lRoming et al.ll2005[ ) began settled ob- 
servations of the field ~80s after the BAT trigger. A 
X-ray afterglow was promptly localized at a position of 
R.A.=05'^55'"2.01^ Dec.=-15°15'19.1'^ (J200 0) with an 
uncertainty of 1.4" (jBeardmore et all I2010I ). while no 
counterpart was observed in the early UVOT unfiltered 
exposures down to a limiting magnitude wh >20.5 (3 cr 
confidence level: [Oates fc Gannizzciir2010[ ). 

Swift data we re analyzed in a standard fashion; we re- 
fer the reader to lEvans et"al[ (|2007ll2010l ) for further de- 
tails. As shown in Fig. 1, the early X-ray afterglow (top 
panel) is characterized by a series of bright X-ray flares 
superimposed on a power law decay (oc t~^'^). Each 
flare can be described by a Fast-Rise Exponential De- 
cay (FRED) profile (sohd line) with 0.04 < At/t < 0.2. 
In the same time interval, the BAT light curve (bot- 
tom panel) shows a long-lasting emission extending up 
to ^^800 s, with several peaks visible in coincidence with 
the X-ray flares (vertical dot-dashed lines). 
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Figure 1. Top Panel: Early XRT light curve of GRB 100728A. 
Bottom Panel: BAT mask-weighted light curve during the X-ray 
flaring activity. Several peaks are visible in correspondance of the 
X-ray flares (vertical dot-dashed lines). 



The postfiare X-ray afterglow decays as a power 
law with slope q;2=1.07±0.05, which steepens to 
a3=1.63±0.07 at t~10 ks. No significant spectral evo- 
lution is observed. The time-averaged photon index is 
r = -2.07 ±0.09. 

By combining the simultaneous BAT and XRT obser- 
vations, we performed a joint spectral analysis of the 
X-ray fiares. We modeled the absorption with two dif- 
ferent components: the former was fixed at the Galac- 
tic value of Nh^IO^i cm^^ (^Kalbcrla ct_al, 2005); the 
latter, representing the absorption local to the burst, 
was fixed to the value of Nh=2.6x10^^ cm~^ derived 
from the late-time (10''-10^s) afterglow spectrum. This 
constraint prevents artificial Njj variations caused by 
the intrinsic spectral evol ution, commonly observe d in 
the brightest X-ray fiares (jButler fc Kocevskil [20071 ) . A 
strong spectral evolution is observed during the first 
100 s, showing a peak energy that softens from 95±15 keV 
during the first flare (from To+167s to To+192s ) to 
less than lOkeV in the following fiares. Excluding the 
first harder episode, the time-averaged spectrum, from 
Tn+ 254s to Tn-|-854s. is well described by a Band func- 
tion ([Band et al.lll993l ) (x^=614 for 477 de grees of free- 
dom, d.o.f.) with a=-1.06±0.11, ^=-2.24±0.02 and a 
peak energy Epk=1.0to;4 keV. 

2.2. Fermi data 

The Fermi /GBM triggered and located GRB 100728A 
at 02:17:31 UT, 53.6 s before the Swift/BAT trigger (see 
§ 12. ip . The GBM light curve shows a complex, multi- 
peaked structure with a duration Tan ^163 s in the 50- 
300 keV energy range (|Kienlin|[20la ). A set of strong 
peaks is visible at To-|-170s, corresponding to the first 
fiares detected by the XRT. No significant emission above 
the highly- variable background level is detected on longer 
timescales. 

The time-averaged spectrum during the Tgo interval, 
from To+15s to To-|-178s can be described with a Band 
function with the following parameters: a=-0.58±0.03, 
P=-2.73to ll and £;pk=264±ll keV (Castor C-statistics 



865 for 351 d.o.f.). A power law function with an ex- 
ponential high-energy cutoff also provides an adequate 
description (C-statistics 885 for 352 d.o.f.). The event 
fiuence (10-1000 keV) in the selected time interval is 
1.181±0.010xl0-'* ergcm-2. The high fiuence of this 
burst generated an Autonomous Repoint Request, which 
caused the Fermi satellite to slew to the GRB position. 

2.2.1. LAT Observations 



An unbinned likelihood analysis (|Abdo et al.l l2009dD 
was used to search the LAT data for emission from 
GRB 100728A. As this study is part of a systematic 
search for high-energy (HE) emission from X-ray fiares, 
a trials factor of 28 for the number of flares considered 
has to be taken into account in evaluating the detection 
significances. 

Depending on the time window of interest the search 
was performed on transient-class data, optimally suited 
for short duration (tens of seconds) signal-limited stud- 
ies, or diffuse-class data, best suit ed for detectin g faint 
emission over longer timescales (jAtwood et al.l [2009) . 
The analysis included LAT events reconstructed within 
15° around the XRT localization (§ 12. ip with energies 
in the 100MeV-50GeV range. The GRB spectrum was 
modeled using a power law. No point source in the 
vicinity of the GRB (within 15°) was bright enough to 
merit inclusion in the background model. The cosmic- 
ray background and the extragalactic gamma-ray back- 
ground were estimated following the method described in 
lAbdo et"al] (|2009d) for the transient-class searches and 
modeled as a single isotropic power law for the diffuse- 
class searches. The Galactic diffuse gamma-ray back- 
ground was described by using the publicly available tem- 
plate produced by the LAT collaboratiori^. The back- 
ground contribution from the Earth's albedo was negli- 
gible since the GRB position was far from the Earth's 
limb during all the time intervals analyzed. Tests per- 
formed with different background models do not show 
any significant change in our results. 

The results of our analysis are summarized in Ta- 
ble [T] A time-resolved search performed on each in- 
dividual fiare did not find any significant excess in the 
LAT data, though a marginal evidence of emission (test 
statistic TS>10, single trial) is present during two of 
the fiares (#3 and #6). A time-integrated search per- 
formed over the whole fiaring interval led to a signifi- 
cant detection (TS=32 for transient-class events, TS=42 
for diffuse-class events). In this time interval the total 
number of transient (diffuse) class events is 191 (29); ac- 
cording to the likelihood analysis the number of events 
associated with the GRB is ^10 (6). The highest en- 
ergy diffuse-class event detected during the fiaring in- 
terval (at ro-|-709s) and in spatial coincidence with the 
source has an energy of 1.68 GeV. The probability of 
the LAT background producing an event with at least 
that energy and during the same interval is Ri7xl0~'*. 
Events of higher energies, tens of GeV, are detected in 
the transient class dataset, but the high background rate 
does not allow us to significantly associate them with the 
GRB. Our best localization of the LAT emission, derived 
from transient-class data analysis, is: R.A.=05'^55™49*', 
Dec.=15°03'18", with a statistical uncertainty of 0.1° 

http; / /fermi. gsfc.nasa.gov / ssc / data/acccss/lat/BackgroundModels.html 
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Table 1 

LAT-Analysis Results 



Time interval 



Test Statistic 



Flux'' 
(lO^S ph cm- 



flare 1 


167-192 





<28 




flare 2 


254-304 





<12 




flare 3 


309-354 


11 


<30 




flare 4 


359-414 





<12 




flare 5 


439-474 





<15 




flare 6 


504-544 


11 


<30 




flare 7 


577-694 





<7 




flare 8 


724-854 





<4 




'i'ime-integratcd ycarch 


pre-flares (prompt) 

post-flares'^ 

X-ray flares 

X-ray flares 


0-167 
854-1654 
167-854 
254-854 


5 

10 
32 
27 


<18 
0.7±0.5 
2.4±1 
2.0±1 


-1.4±0.4 
-1.4±0.2 
-1.3±0.3 


flares 3&6 




22 


9.6±5 


-1.2±0.3 


flaring interval (excluding 3 and 6) 




17 


1.6±1 


-1.6±0.4 



''Fluxes in the 100 MeV - 50 GeV energy band. The quoted errors are at the 68% confidence level. 
Upper limits are at the 95% confidence level and were calculated using the best fit photon index 
rLAT=-l-4. 

''Here the photon index Flat is defined such that dN/dAdEdt oc SI'lat. 
'^From diffuse-class LAT data. 



70 confidence level) and a systematic error of 0.2°. 
This position is consistent with the Swift localization 

(§1211). 

In order to determine whether the LAT emission is 
temporally extended or mainly originated during the 
higher-significance flares (#3 and #6), we performed two 
stacked searches on the transient-class dataset: one ag- 
gregating the data during these two flares and one dur- 
ing the whole flaring period excluding the two flares. 
Emission at a comparable level and with consistent spec- 
tral properties is present during both time intervals (see 
Table [T]), therefore we conclude that the LAT emis- 
sion extends over the whole flaring period. A cross- 
correlation analysis between the LAT (diffuse-class) and 
XRT lightcurves does not detect any significant temporal 
correlation or anti-correlation between the two datasets. 
Similar results are obtained from the analysis of the 
transient-class events. 

As shown in Table [TJ no emission is detected during 
the GRB prompt phase. The resulting upper limit is 
consistent with the extrapolation of the Band spectrum 
to the LAT energy range. Marginal evidence of emission 
(TS«10 for diffuse-class events) is present after the end 
of the observed X-ray flaring activity. 

3. DISCUSSION 

Below we summarize the results that are relevant to 
address the origin of the GeV emission. 

• Significant GeV emission is found in the same in- 
terval where the X-ray fiaring activity is enhanced. 
However, the backgrounds and limited statistics in 
the LAT data do not allow us to search for a one-to- 
one correlation between the GeV emission and the 
single flare episodes. There is marginal evidence of 
GeV emission after the end of the flaring period. 

• The GeV flux is consistent with the extrapola- 
tion of the power law describing the flare spectrum 
above 1 keV. Assuming that an afterglow compo- 
nent is present below the flares and that it has the 



same spectrum observed at later times, it is found 
that the GeV flux is also consistent with the ex- 
trapolation of this putative component. The LAT 
data exhibit a harder spectrum than observed in 
X-rays (see Fig. [5]), though marginally consistent 
(within 3 a) with the X-ray spectral slope. 

The last result suggests that the HE emission can sim- 
ply represent the high energy tail of the synchrotron com- 
ponent. The presence of an additional Inverse Compton 
(IC) component dominating over the synchrotron just 
above ~1 GeV cannot be excluded, and it would be con- 
sistent with the observed flatter GeV spectrum. These 
deductions apply to whichever is the source of electron 
acceleration, internal or external shocks. 




1000 10" 10= 
Energy [keV] 

Figure 2. Spectral energy distribution of the X-ray flares in- 
cluding data from Swift/BKI , XRT and Fermi/LAT. Error bars 
are at 1 cr confidence level. The solid line shows t he be st fit model 
of the joint XRT/BAT spectral fit described in ^ 12.11 and extrap- 
olated to the Fermi/LAT energy range. As discussed in the text, 
the Fermi/hNT detection is consistent within 3 ct with the extrap- 
olation of the model. 
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The discovery of HE emission in a time frame of vig- 
orous flaring activity in X-rays lead us to consider firstly 
the association of the GeV emission with X-ray flares. 
We now discuss this scenario. Given the large num- 
ber of flares, we can exclude that a delayed external 
shock is the domina nt process originating the X-ray flares 
(jGalh fc Piroll2Q07[ ). In fact, in this model only a single 
outstanding flare, corresponding to the onset of the af- 
terglow from a long duration central engine, is produced. 
We thus consider internal shocks from a lon g-lasting rel- 
ativis tic outflow as the source of flares (e.g. I Zhang et al.l 
I2OOI . 

In order to allow the GeV emission to be observed we 
require two conditions. First, the source has to be op- 
tically thin for pair production. By computi ng the op- 
tical depth for e+e" production ([Lithwick fc Sari 2001 
by photons of energy EceV in GeV on the X-ray to GeV 
power law component observed at about 300 s we derive 
a lower limit on the Lorentz factor; 



r > r K,'\V\ K^/® r)i/3 



1/3 



(1) 



where we specialized the equation for a photon index 2. 
The tightest constraints on F are derived from the short- 
est time scale for variability i„ that can be associated 
to the relativistic flow produced by the central source. 
In the scenario of late internal shock X-ray flares and 
the GRB prompt emission are both related to the cen- 
tral engine activity. We therefore consider that a vari- 
ability timescale of ms typical of the prompt phase is 
a reasonable possibility. In this case a single flare is 
produced by the superposition of several internal shock 
events from a relativistic wind active for the whole dura- 
tion of the flare. By assuming t„ = 10^"^ s, a time scale 
similar to that characterizing the prompt phase, one de- 



rives F. 



77 



115 Eq^y for a typical redshift z = I. On 
the contrary, if the flare is associated to a single internal 
shock event, i.e. the interaction of two shells, then ty 
should be equal to the flare duration, i.e. about 100 s. In 



this case F, 



20 E, 



1/6 



77 ^"-^Gey 

The lower boundary on the Lorentz factor derived 
above for the flaring phase encom passes the range of val- 
ues typical of the prompt phase ([Lithwick fc Sarill2001l : 
ILiang et al.ll2QlCil ) , consistent with the notion that a long 
duration relativistic outflow with a Lorentz factor of the 
order of wlOO is producing both the prompt emission 
and the flares. In principle the parameters describing 
the relativistic shocks (eg, es, L, F, t^) can all be time de- 
pendent, i.e. be different during the prompt and late 
flaring phases. On the other hand, one wishes to re- 
duce the number of variable parameters (Occam's razor). 
It goes beyond the scope of this paper to find the best 
self-consistent internal shock model reproducing both the 
prompt and X-ray flaring phases. We just note the fol- 
lowing. The model should be able to reproduce a peak 
energy that shifts from the «100 keV region during the 
prompt phase to the keV range observed during X-ray 
flares. Recalling that the peak of the synchrot ron spec- 
trum is given by fe.g. iZhang fc Meszarosl[200^ : 



it follows that the decrease of the luminosity L from the 
prompt to the flare phase already accounts for a decrease 
of the peak energy by a factor of 20, with all the other 
parameters remaining constant. The further reduction 
that is needed can be obtained e.g. by the very reason- 
able assumption that the magnetic field weakens at the 
larger radii where flares are produced or by a smaller 
contrast in the Lore ntz factor between colliding shells 
(iBarraud et al.|[2005l ). 

The second condition is derived by requiring that the 
maximum energy at which the electrons are accelerated 
is large enough to produce photons of energy E via syn- 
chrotron radiation: 



F > 60 



l + z 



il+Y)EGeV, 



(3) 



where Y is the Compton y parameter. This equation 
gives a condition on F comparable to that derived from 
Eq. [TJ In conclusion we find that both the prompt emis- 
sion and the later X-ray flares and HE emission can be 
explained by internal shocks produced by a long dura- 
tion central engine with a Lorentz factor of w 100 and 
decreasing luminosity. 

This simple internal shock model predicts an emission 
that is co-spatial and simultaneous in the X-ray and 
GeV ranges. On the other hand, we find a marginal 
evidence of delayed HE emission. This is naturally pre- 
dicted when the X-ray photons, produced by internal 
shocks at smaller radii, are upscattered to GeV energies 
via IC by t he electron population of the forward shock 
(|Wang et a l. 2006; Fan et al. 2008). 

Finally, given the quality of the present dataset, we 
cannot exclude the possibility that the GeV emission 
is actually related to an afterglow underlying the X-ray 
flares. This requires that the afterglow onset takes place 
before 200 s. Such condition is satisfied when the Lorentz 
factor of the relativistic fiow at the beginning of the de- 
celeration phase is: 



Fea;* > 



1/4 



ISM, 
Wind, 



(4) 



3/2 l/2.1/2p-2,-l 



(2) 



where A* is the density scaling factor in units 
of 5x10" gcm"^ In other GRBs, the HE emis- 
sion has been indeed associated to the f o rward 
shock synchrotron emission (Ghir landa eTall IMl 
IKumar fc Barniol Duranll2010l) , though the shorter time 
scales observed require much larger values of Text than 
those derived here. We further explore the external shock 
scenario in GRB 100728A by analyzing the late time X- 
ray behavior. The afterglow spectral and temporal laws 
are bound to obey specific r elations (the so- called clo- 
sure relations, e.g. IZhang &: Meszaroil2004D . that de- 
pend upon the density profile of the external medium, 
the jet opening angle, and the relative position of the 
typical frequencies of the synchrotron spectrum with re- 
spect to the observed range. Within the simple external 
shock model, the closest solution envisages a jet with a 
rather narrow opening angle of ~l-2deg expanding in 
a medium with a wind-like density profile, though the 
lack of multi-wavelength afterglow observations does not 
allow us to firmly characterize the circumburst environ- 
ment. This scenario is consistent with the lack of spectral 
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variations before and after the break at 10 ks, albeit in 
a wind-like medium a jet transitio n is expected to take 
place on much longer time scales ()Kumar fc Panaitescul 
[2000l ). In this scenario the cooling frequency falls below 
the X-ray band, and the early GeV emission (if after- 
glow) likely belongs to the same synchrotron regime. In 
this case, the LAT emission should display a similar de- 
cay slope of ~1.07 and a photon index of ^-2.07, softer 
than the observed value of -1.4±0.2 (la) but still con- 
sistent within the large uncertainty. 

4. CONCLUSION 

GRB 100728A is the second case to date with si- 
multaneous Swift and Fermi observations. High-energy 
gamma-rays are detected by the Fermi/h AT until 850 s 
(TS=42) and possibly continuing until 1600s (TS«10). 
Very interestingly the early X-ray afterglow exhibits an 
intense and long-lasting flaring activity, visible both in 
BAT and XRT. Although an afterglow origin of the GeV 
emission cannot be excluded, the presence of bright X-ray 
flares unveiled by Swift observations opens the possibility 
that a prolonged central engine activity is powering the 
temporally extended HE emission observed in this burst. 

Within the internal shock scenario a relativistic out- 
flow with a Lorentz factor of w 100 and decreasing lumi- 
nosity can explain the prompt emission, the later X-ray 
flares and HE emission. The presence of a delayed HE 
emission naturally arises from IC scattering of low-energy 
flare photons off the relativistic electrons at the external 
forward shock radius. 
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